APS/123-QED 



Equality of Bond Percolation Critical Exponents for Pairs of Dual Lattices 



o 
o 

(N 



Matthew R. A. Sedloclfl and John C. Wiermai|!] 

Applied Mathematics and Statistics Department, The Johns Hopkins University, Baltimore, Maryland 21218 

For a certain class of two-dimensional lattices, lattice-dual pairs are shown to have the same bond 
percolation critical exponents. A computational proof is given for the martini lattice and its dual 
to illustrate the method. The result is generalized to a class of lattices that allows the equality 
of bond percolation critical exponents for lattice-dual pairs to be concluded without performing 
the computations. The proof uses the substitution method, which involves stochastic ordering of 
probability measures on partially ordered sets. As a consequence, there is an infinite collection 
of infinite sets of two-dimensional lattices, such that all lattices in a set have the same critical 
exponents. 
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I. INTRODUCTION 

The behavior of a percolation system near the critical 
probability is often expressed in terms of critical expo- 
nents. Their values are believed to depend only on the 
dimension of the lattice, rather than the structure of the 
lattice itself. This is known as the universality hypothe- 
sis. Using scaling and hyperscaling relations, values for 
the critical exponents in two dimensions have been pro- 

Sised and supported by substantial simulation evidence 

The works of Kesten [4I and Wierman [3| give mathe- 
matical progress toward showing the universality hypoth- 
esis in two dimensions; Kesten was able to show relation- 
ships among critical exponents assuming that the lim- 
its defining them exist, while Wierman was able to es- 
tablish the equality of certain bond percolation critical 
exponents for two pairs of dual lattices: the triangular- 
hexagonal pair and the bowtie-dual pair. 

Very little has been proved about critical exponents in 
two dimensions. Kesten Q shows that specific functions 
can be bounded above and below by powers of |p — Pc| 
for lattices in which the "horizontal and vertical direction 
play symmetric roles." In particular, these results do not 
rely on the calculation of the percolation threshold, Pc, of 
a given lattice, but rely instead on other properties of the 
lattice. Wierman's result for the triangular-hexagonal 
pair and the bowtie-dual pair rely heavily on the ability 
to calculate the percolation threshold for these lattices 
using the star- triangle transformation. As a result, little 
use has been made of this method since the exact per- 
colation threshold is known for relatively few lattices. A 
notable advance by Smirnov and Werner Q is the deter- 
mination of the values of the critical exponents of the site 
percolation model on the triangular lattice. 

The recent work of Ziff and ScuUard [1,0, 11] generalizes 
the star-triangle transformation to find the exact perco- 



lation thresholds for the martini lattice and its variants, 
the A lattice and the B lattice. Using this idea, Wier- 
man and Ziff Q have since identified an infinite class 
of lattices-those constructed from an infinite connected 
planar periodic 3-uniform hypergraph with one axis of 
symmetry, using a generator which is a finite connected 
planar graph with three boundary vertices-for which the 
exact bond percolation threshold can be determined. By 
applying Wierman's method to the martini lattice and its 
dual, this paper gives a third lattice-dual pair with the 
same bond percolation critical exponents. More impor- 
tantly, the method can be extended to the class identified 
by Wierman and Ziff, so that infinitely many lattice-dual 
pairs are shown to have the same values for their critical 
exponents. Although the ability to calculate the perco- 
lation threshold is essential for proving the result for lat- 
tices in the class mentioned, the result does not depend 
on the percolation threshold of a lattice. 

Section II provides some background information, de- 
scribes the class of lattices, and explains the method. 
Section III uses the method of Wierman Q to show equal- 
ity of the critical exponents for the martini lattice and its 
dual. Section IV shows that the method of Wierman can 
be extended to lattices constructed from an infinite con- 
nected planar periodic 3-uniform hypergraph with one 
axis of symmetry, using a generator which is a finite con- 
nected planar graph with three boundary vertices. Thus, 
without actually going through the calculations, we can 
conclude that a lattice belonging to this class and its dual 
have the same bond percolation critical exponents. Sec- 
tion V shows that the result includes other lattices that 
can be transformed into lattices in the class and discusses 
the implications for the equivalent site problems. 



II. PRELIMINARIES 



A. Percolation Functions and Critical Exponents 
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In the bond percolation model on an underlying graph 
G, bonds between vertices are labeled open with proba- 
bility p, < p < 1, independent of all other bonds. A 
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bond that is not open is said to be closed. The probabil- 
ity measure and expectation operator corresponding to a 
given p are Pp and Ep respectively. Several interesting 
quantities arise in this model. For any given vertex v in 
G, let Cv be the set of bonds in the connected component 
containing v. Then, | Cv \ gives the number of bonds in the 
connected component containing v. The main quantities 
considered in this paper are the percolation probability 
functions, Oy{p) = Pp[|Ci,| = +00], and the mean finite 
cluster size, Xv{p) = Ep[\Cv\;\Cv\ < cxd]. The two-point 
connectivity function and the correlation length are also 
of interest and are defined in section V. 

It is believed, though not proved, that most quantities 
of interest in percolation theory behave as powers of \p — 
Pc\ as p approaches Pc- These powers are called critical 
exponents. The notation 

is used to mean that 

lim . r = C- 

P^Pa log \p~Pc\ 

It is not known that these limits exist, so we use the 
superscripts -I- and — on the exponent to denote limsup 
and liminf respectively. The power laws we consider are 

0{p) w {p-pcf, for p>pc, 

for some < /3 < 1, and 

Xfip) ~ b-Pcr'^ 

for some 7 > 0. The power laws for the two-point connec- 
tivity function, the correlation length, and other quanti- 
ties are given in section V. 

B. The Class of Lattices 

We now describe the class of lattices, introduced by 
Wierman and Ziff Q, for which our results are valid. We 
will refer to it as the "martini class" of lattices. The lat- 
tices are constructed by placing copies of a planar graph 
called a generator into a self-dual 3-uniform hypergraph 
arrangement. The concepts of generator, 3-uniform hy- 
pergraph, and self-duality in this context are explained in 
the following subsections. 

1. Self-dual Hypergraph Arrangements 

The concept of hypergraph is a generalization of the 
concept of graph, in which hyperedges connect sets of 
vertices rather than edges connecting pairs of vertices. 
Given a set V of vertices, a hyperedge H is a subset of 
V. A hyperedge H is said to be incident to each of its 
vertices. A hyperedge containing exactly k vertices is 
called a /c-hyperedge. 



A hypergraph is a vertex set V together with a set of 
hyperdeges of vertices in V. A hypergraph containing 
only fc-hyperedges is a fc-uniform hypergraph. A hyper- 
graph is planar if it can be embedded in the plane with 
each hyperedge represented by a bounded region enclosed 
by a simple closed curve with its vertices on the bound- 
ary, such that the intersection of two hyperedges is a set 
of vertices. 

In this article, we will consider only 3-uniform hyper- 
graphs. For convenient visualization, we will represent 
each 3-hyperedge in the plane as a shaded triangular re- 
gion bounded by a slightly concave triangular boundary. 
This allows us to neglect the detailed structure of our 
generators when we consider arranging them in the plane 
to form a connected structure. 

In order to construct exactly-solvable lattice graphs for 
bond percolation models, we will consider infinite con- 
nected planar periodic 3-uniform hypergraphs. A planar 
hypergraph H is periodic if there exists an embedding 
with a pair of basis vectors u and v such that H is invari- 
ant under translation by any integer linear combination 
of u and v, and such that every compact set of the plane 
is intersected by only finitely many hyperedges. 

If a hypergraph H is planar, we may construct its dual 
hypergraph H* as follows. Place a vertex of H* in each 
face of H. For each hyperedge e of H, construct a hy- 
peredge e* of H* consisting of the vertices in the faces 
surrounding e. Note that if the hyperedge of is a 3- 
hyperedge represented by a triangular region, and each of 
its boundary vertices is in at least two hyperedges, then 
the dual hyperedge is a 3-hyperedge also, represented by 
a "reversed triangle." 

Two hypergraphs are isomorphic if there is a one-to- 
one correspondence between their vertex sets which pre- 
serves all hyperedges. A hypergraph is self-dual if it is 
isomorphic to its dual. If, in addition, the hypergraph 
is 3-uniform, this corresponds to the term triangle-dual 
used by Ziff and ScuUard. To illustrate. Figure 1 pro- 
vides two examples of infinite connected planar periodic 
self-dual 3-uniform hypergraphs mentioned by Ziff and 
ScuUard. 

Wierman and Ziff 9] note that the reversed triangles 
must be connected in a specific manner to create a self- 
dual arrangement, and provide an example to illustrate 
that simply reversing each triangle is not sufficient. 



2. Generators and Duality 

A generator is a finite connected planar graph embed- 
ded in the plane so that three vertices on the infinite face 
are designated as boundary vertices A, B, and C. (For 
future research, the term generator may be defined more 
generally. However, in this paper, we restrict attention 
to generators with three boundary vertices.) 

Given a generator G, we construct its dual generator 
G* by placing a vertex in each bounded face of G, and 
three vertices A*, B* , and C* of G* in the infinite face 



3 




FIG. 1: Self-dual hypergraph arrangements illustrated in Ziff 
and ScuUard. We refer to the upper drawing as the triangular 
arrangement and refer to the lower drawing as the bow-tie 
arrangement. 



of G, as follows: The boundary of the infinite face can 
be decomposed into three (possibly intersecting) paths, 
from A to B, B to C, and C to A. The infinite face 
may be partitioned into three infinite regions by three 
non-intersecting polygonal lines starting from A, B, and 
C. Place A* in the region containing the boundary path 
connecting B and C, B* in the region containing the 
boundary path connecting A and C, and C* in the region 
containing the boundary path connecting A and B. A*, 
B* , and C* are the boundary vertices of G* . 

For each edge e of G, construct an edge e* of G* which 
crosses e and connects the vertices in the faces on oppo- 
site sides of e. If e is on the boundary of the infinite face, 
connect it to A* if e is on the boundary path between B 
and C, to B* if e is between A and C, and connect it 
to C* if e is between A and B. (Note that it is possible 
for e* to connect more than one of A* , B* , and C*, for 
example, if there is a single edge incident to A in G, its 
dual edge connects B* and G* .) 

Note that G* is not the dual graph of G, which would 
have only one vertex in the infinite face. The three ver- 
tices A*, B* , and G* will correspond to separate faces of 
the lattice generated from G. 



3. Constructing a Dual Pair of Lattices 

Given a planar generator G and a connected periodic 
self-dual 3-uniform hypergraph 7i, a dual pair of peri- 
odic lattices may be constructed as follows: Construct 
a lattice graph Lg,h by replacing each hyperedge of H 
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A* 

FIG. 2: Solid lines represent a 3-hyperedge with boundary 
vertices A, B, and C. Dashed lines represent the "reversed" 
or dual hyperedge, with its boundary vertices A* , B* , and C* 
labeled in the proper positions. We say that A* is opposite 
A, B* is opposite B, and C* is opposite C. 

by a copy of the generator G, with the boundary ver- 
tices of the generator corresponding to the vertices of 
the hyperedge, in such a manner that the resulting lat- 
tice is periodic. This is always possible, by choosing the 
orientations of the generator in one period of the hyper- 
graph, and extending the choice periodically. (However, 
for a generator without sufficient symmetry, it may be 
possible to replace hyperedges in a way that produces a 
non-periodic lattice, so some care is needed.) 

We now construct a lattice Lew as follows: Con- 
struct the embedding of the dual hypergraph 7i* in the 
plane, in which every hyperedge of 7i is reversed. Replace 
each hyperedge of Ti* by a copy of the dual generator G* , 
embedded so that it is consistent with the embedding of 
G, that is, in all hyperedges boundary vertex A* in G* 
is opposite vertex A in G, B* is opposite B, and G* is 
opposite G, and each edge of G* crosses the appropriate 
edge of G. (See Figure 2.) This results in a simultane- 
ous embedding of Lg*,w and Lcn- example of the 
construction for a particular generator is illustrated in 
Figure 3. 

The constructions of the two lattices both produce a 
planar representation of the resulting lattice. From the 
simultaneous embeddings of the two lattices, it is seen 
that Lg'^h' is the dual lattice of Lq^-h^ since there is a 
one-to-one correspondence between vertices of one and 
faces of the other, and a one-to-one correspondence be- 
tween edges, which are paired by crossing. 



C. Substitution Method 

Wierman introduced the substitution method to 
find bounds on bond percolation critical probabilities for 
certain lattices. Its application to the martini class of 
lattices is described here. Let Q he a, lattice constructed 
by placing copies of a generator in a self-dual 3-uniform 
hypergraph arrangement. A boundary vertex is a vertex 
that is in more than one copy of the generator. Let G be 
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FIG. 3: The construction of lattices based on a specific gen- 
erator. Top: The generator, the duality relationship, and the 
dual generator. Middle: The lattices based on the generator 
and the triangular hypergraph arrangement. Bottom: The 
lattices based on the generator and the bow-tie hypergraph 
arrangement. 



the generator and suppose that its 3 boundary vertices 
are A, B, and C. Any configuration of open and closed 
bonds on G gives a partition of the boundary vertices 
into connected open clusters of vertices. Vertical bars 
between the boundary vertices are used to denote such 
a partition. For example, AB\C denotes the partition 
where A and B are in the same open cluster and C is in 
a separate open cluster. 



1. Partially Ordered Sets 

A partially ordered set, or briefly, a poset, consists of a 
pair (5, <), where 5 is a set and < is a binary relation 
with the following properties: (1) < is reflexive, i.e., s < s 
for all s € S. (2) < is antisymmetric, i.e., for all s and t 
in 5, if s < i and t < s, then s = t. (3) < is transitive, 
i.e., for all s,t, and u in S, ii s <t and t < u, then s < u. 

In a partially ordered set {S, <), two elements of S, s 
and t, arc comparable ii s < t or t < s, and are incompa- 
rable otherwise. A partially ordered set is a total order 
if every pair of its elements are comparable. For exam- 
ple, the standard relation < on a set of real numbers is a 
total order. Partially ordered sets generalize the concept 
of total order by allowing incomparability. A common 
example of a partially ordered set that is not totally or- 
dered is the set of all subsets, or power set, of a set A of 
two or more elements with the relation of set inclusion, 
since some pairs of subsets are incomparable. 

It is useful to have a visual representation of a partially 




FIG. 4: A Hasse diagram of the partially ordered set of par- 
titions of three boundary vertices, ordered by refinement. 



ordered set. To efficiently represent a poset, it is useful 
to define the cover relationship. For s,t £ S, s is covered 
by t ii s < t and there is no element u & S, unequal to 
both s and t, such that s < u < t. The Hasse diagram 
of the partially ordered set {S, <) is obtained by placing 
a point in the plane for each element of S, taking care 
to put the point for ,s below the point for t whenever 
s < t, and connecting the points for s and i by a line 
segment if and only if s is covered by t. Note that all 
other comparability relationships can be obtained from 
the covering relationships by following monotone paths in 
the Hasse diagram, so it is not necessary to complicate 
the diagram by including line segments for all pairs of 
comparable elements. 

The set of boundary partitions of a generator is a 
partially ordered set, with, for boundary partitions tt 
and a, TT < a whenever every cluster of a can be de- 
composed into clusters of tt. In the previous example, 
^|i?|C < AB\C, since A\B is a decomposition of AB 
and C is a trivial decomposition of C. Equivalently, ev- 
ery cluster of tt is entirely contained in a single cluster of 
a. Another equivalent condition is that every cluster of 
a is a disjoint union of clusters of tt. 

Figure 4 illustrates the Hasse diagram of a partially or- 
dered set with 5 elements which is used throughout this 
article. The three middle partitions are pairwise incom- 
parable, while all other pairs of partitions are compara- 
ble. 



2. Stochastic Ordering 

Stochastic ordering is used to compare two probability 
measures defined on the same partially ordered set. We 
will consider probability measures which are defined on 
the set of boundary partitions of a generator in terms 
of a bond percolation model on the generator. Consider 
a generator G and a bond percolation model on G with 
parameter p. A configuration is a designation of each 
edge of G as open or closed. In the percolation model, 
a configuration with k open edges and I closed edges has 
probability — p)'. A boundary partition is a union 
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of configurations, with its probability being the sum of 
the probabihties of its configurations. Let denote the 
probabihty measure on boundary partitions of G gener- 
ated by the bond percolation model with parameter p. 

Let 5 be a partially ordered set and let U C S. U 
is called an upset if for all / and g in S, f < g and 
f & U imply that g G U. UP and Q are two probability 
measures on S, we say P is stochastically smaller than 
Q, denoted P <s Q, if P{U) < Q{U) for all upsets U of 
S. This concept of stochastic ordering is the appropriate 
comparison of two probability measures on a partially 
ordered set when applying the substitution method Q . 

In addition to a lattice Q, suppose another lattice Ti. 
is constructed by placing copies of a generator H in 
the same self-dual 3-uniform hypergraph structure as G 
The substitution method studies the effect on connec- 
tion probabilities of replacing the generator H oi H by 
the generator G in Q. 

Preston [Tl| established the equivalence of stochastic 
ordering and coupling for probability measures on finite 
partially ordered sets. For simplicity here, we will spe- 
cialize to our percolation setting. If P^ <s P^ for two 
generators G and H, then there exist percolation models 
on G and H with parameters p and q respectively which 
are dependent on each other in such a manner that when- 
ever boundary vertices are connected by open edges in G 
they are also connected by open edges in H . (Note that 
each percolation model corresponds to independently 
open edges in its generator, so is determined by its lat- 
tice connectivity, but that the realizations of the two per- 
colation models are stochastically dependent upon each 
other.) In particular, the stochastic ordering P^ <s Pq 

implies that 9^ {p) < e^{q) and x^{p) < X^{q)- We use 
these consequences in section III, where solving for (5(e) in 
the equations Pp^^^{U) = P^^g(^^-^{U) for all upsets U al- 
lows us to obtain a stochastic ordering on the probability 
measures for a lattice and its dual near their percolation 
thresholds. 



III. EXAMPLE: THE MARTINI LATTICE AND 
ITS DUAL 

In this section, we show that the critical exponents of 
the bond percolation models on the martini lattice and its 
dual, referred to as the K4 lattice, are equal. The gener- 
ators for the martini and lattices are shown in Figure 
5. The method used is that introduced by Wierman Q to 
show the equality of the bond percolation critical expo- 
nents for the triangular and hexagonal lattices and also 
for the bowtie lattice and its dual. There are no other 
cases where Wierman's method has been used to show the 
equality of bond percolation critical exponents for dual 
lattices. This example is used as a reference in section IV, 
where the main result is shown. Precisely, we show that 
P+{M) = p+{K), 13- {M) = p-{K), -f+{M) = -f+{K), 
and 7 (M) = 7 {K), where M and K refer to the mar- 
tini and Ki lattices respectively. So, if the limit defining 



the critical exponent /? exists for either lattice, it exists 
for both lattices and /3(M) = (3(K). Similarly, if the 
limit defining the critical exponent 7 exists for either lat- 
tice, it exists for both lattices and 7(M) = "i{K). As in 
Wierman's paper, the substitution method is used in this 
proof. 




FIG. 5: Generators for the Martini and K4, Lattices, on the 
left and right, respectively. 

Ziff [gj was able to find the exact percolation thresh- 
old for the martini and A'4 lattices. A derivation of this 
threshold is summarized here. Let each bond in the K^^ 
lattice be open with probability p. Then, the partition 
probabilites for the K/^ lattice can be calculated by con- 
ditioning on the bonds in the interior star: 

P^[ABC] = 3p2 + 5p3 _ 18/ + i5p5 _ 

P^[AB\G] = P^[AC\B]^ P^[A\BC] 

= p-p^ -5p^ + Up^ - + 
P^[A\B\G] ^l-3p+ lOp^ - 15/ + 9p5 - 2/ 




FIG. 6: The K4, Lattice 

Similarly, letting each bond be open in the martini 
lattice with probability q, the partition probabilities for 
the martini lattice can be calculated by conditioning on 
the bonds in the triangle: 

P^^ABC] = 3q^ - 2q^' 

P^^[AB\C] = P^^[AG\B]^P^^[A\BC] 

= q^ + q^ ~ Aq^ + 2q^ 
P^^[A\B\G] = 1 - 3g3 - 3g4 + 9g5 _ 4^6 
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FIG. 7: The Martini Lattice 



Using the substitution method, we replace each mar- 
tini generator by the K4 generator and set P^[ABC] ~ 
P,^[ABC], Pp^[AB\C] - P,^nAB\C], and Pp^[^|B|C] = 
Pq*'^[j4|i?|C]. The second equation is satisfied when p — 
1 — q, which is a consequence of duahty of the two lattices. 
The first and last equations are redundant, and give the 
following: 

{2q^ -l){q^ -3q^ + 2q^ + 1) = 

which has roots a,t q — We thus have that the crit- 

ical probabilities are Pc{M) = and Pc{K) = 1 — 
as shown before by Ziff [6| . 

To show equality of the bond percolation critical ex- 
ponents, we will establish a stochastic ordering between 
the probability measures P^ near p — Pc{K) and P^^ 
near q = qdM). Specifically, we find that 



M 

gc + -93e 



< 



pK ^ 



sP, 



by showing that 



nt + .93e(t/) < P/^+,(t/) < P*Vl.07e(t/) 

for all upsets U of the poset of boundary partitions. This 
allows us to make conclusions about the behavior of the 
percolation probability near criticality. Proceeding, we 
equate the upset probabilities for the martini and 
lattices, keeping in mind that they are equal at criticality 
as a result of duality. The upset probability equations are 
given by: 

3p2 -I- 5p3 „ 18p4 ^ i^p5 „ 4p6_j_ 

k[p -p^ - 5p3 llp4 _ ^ 2/] 
3q5 - 2q^ + k[q^ + q^ - 4(7^ + 2q^] 

for fc = 0, 1, 2, 3, where k denotes the number of middle 
partitions in the upset. In studying critical exponents, we 
consider the behavior of the system near the percolation 
threshold, so we consider perturbations near pc and qc- 
Thus, we want to solve for 5 = 5{e) in the following: 

3(pc + e)' + 5(pc + e)^-18(pc + e)^ + 15(pc + e)5-4(pe + e)6 
+k[{p^ + e) - {pc + - 5{pc + 
+n{pc + e)* - 8{pc + ef + 2{p, + e)^] 



= 3(<z, + 5)^ ~ 2{q, + 5f 
-k[{q, + 6f + {q, + (5)4 - 4(q, + 5f + 2{q, 



5f 



for k = 0,1,2,3. 

Noting that the constant terms cancel out by the upset 
probability equations and letting 

g{p) = (6p + 15p2 - 72p3 + 75p4 _ 24p5) 

+fc(l -2p- Ibp^ + Up^ - 40/ + 12p5) 



and 

h{q) = 15/ 
we see that 



I2q^ + k{3q^ 



iq^ - 20(7" 



o{e). 



12g^ 



.2929, the 



Using the approximation pc = 1 — .7071 
coefficients of e are approximately: 

1.0655 for fc = 
1.0214 for fc = 1 
.9791 for A: = 2 
.9385 for fc = 3 



Thus, for sufficiently small e, .93e is smaller than the 
coefficient in the solution (5(e) for all four equations, and 
1.07e is larger than the coefficient in the solution (5(e) for 
all four equations. Since the upset probability functions 
are increasing functions of the parameters, for sufficiently 
small positive e, we have 

3(pc + e)' + 5(pc + e)^-18(pc + e)" + 15(pc + e)''-4(pc + e)' 



iPc 



+niPc + e)4 - 8{p, + ey + 2{p, + e)^] 

> 3(gc + ■93e)5 - 2{qc + .93e)'^ 
+k[{q^ -h .93e)3 + [qc + Mef 
-A(q^ + mef + 2{qc + m€)% 

for fc = 0, 1, 2, 3, and 

Z{pc + ef + h{pc + €f-lS{p, + e)^ + lb{pc + ff-'^{Pc + ef 



+fc[(pc + 
+ll(p, + ef - 8{pc + ef + 2{p, + 



(Pc 



HPc 



< 3{qc + l.OTef - 2{q, 
+k[{qc + 1.07e)3 + {qc 
-4{qc + 1.07e)5 + 2{qc - 



F 1.07e)6 
- 1.07e)4 
1.07e)S], 



for k = 0,1,2,3. 

Each of these sets of inequalities establishes a stochas- 
tic ordering result. Thus, for sufficiently small e > 0, 



pM ^ pl\ ^ pl\ 

^gc + .93c — -S pc+c — -S gc-l-1.07c 

Using the result from Preston [ill, have that, for 
p > Pc{K) and sufficiently close to Pc{M), 

e'\qc + m{p-pc)) < e'^ip) < e^Hq^ + im{p-pc)) 
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Taking logarithms and dividing by log |p— Pc| throughout, 
we see that 

hmsup. ■ : ^ hmsup — — 

qlqo log|g-gc| Pip, log|.93(p-pc)| 

logg^^(g, + .93(p-pc)) 

= iim sup — 7-— ; ;— 

piPc log(.93) + log|p-pc| 

logg^^(g, + .93(p-pc)) 
= hm sup j j 

Pip. log |P- Pel 

> lim sup j - 

plPa log |P- Pel 

Reversing this argument using q = qc + 1.07(p — Pc) 
gives: 

logg^(p) ^, logO^Hq) 
lim sup j r > lim sup : 



pjPc 



log \p - Pc 



qlqa 



log|g-9c 



Consequently, P^{K) = /3+(M) and, by changing 
limsup to liminf where it appears above, we have that 
p-iK)=p-iM). 

By the same reasoning, we also have that for suffi- 
ciently small e > 0, 

Thus, for p < Pc{K) with p sufficiently close to pc{M), 
X^'iqc - l.OTe) < x^'iPc - e) < X.^' {Qc - .93e). 

Letting e i gives that j+{K) = 7+(Af) and ^^(K) = 
7-(M). 



IV. POSITIVE COEFFICIENTS FOR e IN THE 
5{e) FUNCTIONS 

In this section we generalize the result to lattices 
constructed from an infinite connected planar periodic 
3-uniform hypergraph with one axis of symmetry, using 
a generator which is a finite connected planar graph 
with three boundary vertices. See Wierman and Ziff 
[sl for the construction of such lattices. As in section 
III, the method used to show equality of the critical 
exponents is valid if the coefficients of e in the (5(e) 
functions are positive and finite. In this section we show 
that the coefficients for lattices in the class are positive 
and finite. In what follows, G will denote the generator 
of a lattice in the class of interest and D will denote its 
dual generator. 

Calculate the coefficient of e in the (5(e) function as fol- 
lows. First, equate the upset probabilites for the lattice 
and its dual, setting 



P^[ABC] - 
'[ABC] 



p 



-jP^[AB\C]^ 
^]P^[AB\C] 



kP^[AC\B] 4 
" [AC\B] 



kP^ 



IP^[BC\A] ^ 
[BC\A], 



where j, fc, and I take the values or 1. Note that this is 
more general than the martini and lattices example 
in that symmetry of the generator is not assumed. The 
critical probabilities for the lattice and its dual are then 
calculated using these equations. Adding e to p and 5 to 
q in these equations shows how the two functions behave 
around the critical probabilities. 

Notice that the numerator of the coefficent of e in 
the martini and lattices example is the derivative 
of P^[ABCf\ + kP.^[A\BC] and the denominator 
is the derivative of P^[ABC] + kP^[A\BC]. This 
holds for all generators in the class. That is, the 
numerator of the coefficent of e is the derivative of 
P^[ABC] + ]P^[A\BC] + kP^[AB\C] + IP^[AC\B] 
and the denominator is the derivative of 
Pl^[ABC] + ]Pl^[A\BC] + kP!;'[AB\C] + lPl^[AC\B] 
for any generator in the class. This can be seen by 
the following reasoning: Suppose the generator G has 
n bonds. Then, since the lattice's generator and its 
dual have the same number of bonds, the expressions 
P^[ABC] + jP^[A\BC] + kP^[AB\C] + lP^[AC\B] and 
P°[ABC] + ]P^[A\BC] + kP^[AB\C] + IP!^[AC\B] 
are polynomials in p and q respectively with degree no 
larger than n. So, we can write 

P^[ABC] + ]P^[A\BC] + kP^[AB\C] + IP^[AC\B] 

n 

= akP^ = g{p) 

k=0 

and Pq^[ABC] + jP^ [A\BC] + kP^ [AB\C] + IP^ [AC\B] 

n 

^^bkq'' - h{q) 

Equating these, and adding e and S to p and q respec- 
tively, yield the following: 

n n 

J2bk{q + 5)'' = Y,^k{p + ^f■ 
k=o k=Q 

Using binomial expansions, we can write this as: 



J2bkq'' + Sj2kbkq''-'+o{S) 

n n 

^ akP^ + e ^ kakP^-^ + o(e) 



fe=0 



k=l 



P 

ipF' 



Evaluating this at p = Pc and q — q^ the constant 
terms on each side of this equation cancel as a result of 
the equality of the probability measures for dual lattices 
at criticality. Moreover, the coefficients of 6 and e are 
seen to be the derivatives of their respective upset prob- 
ability functions. We thus have that 
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To show that the coefficient of epsilon is positive and 
finite, it suffices to show that the derivatives of the upset 
probabihty hmctions in both the lattice and its dual are 
positive on the interval (0,1). Since the upset probability 
functions are polynomials, their derivatives exist and are 
finite. If the derivatives are shown to be positive, then 
they will be positive at p,. and the coefficient of e will be 
positive, completing the argument. 

We now prove that the derivatives of the partition 
probability functions for a lattice in the martini class 
and its dual are positive. Let T be a minimal connected 
subgraph of the generator G that contains A, B, and C. 
Let E{T) denote the edge set of T and let m = \E{T)\. 
T is clearly a tree, so there is an unique path from A to 
B, from A to C, and from B to C. More importantly, 
removing any edge from T makes connectivity between 
A, B, and C impossible in T. 

We shall use the following notation for the proof. The 
superscript on the probability measure will indicate the 
structure for which the boundary vertices can be con- 
nected through; G will indicate the generator of the lat- 
tice, T will indicate the tree defined above, and N will 
indicate through G but not through T. (The meaning of 
this will become clearer in the body of the proof.) Where 
no superscript appears, the statement is true for all three 
structures. The subscript will denote the probability of 
an edge being open in the graph. Where two subscripts 
appear, separated by a comma, the first subscript gives 
the probability that ei is open, and the second subscript 
gives the probability that each edge other than ei is open, 
ei may be defined differently in different cases, but ei will 
always be an edge in T. Where no subscript appears, the 
statement is true for any subscript. It will be helpful to 
condition on ei being open or closed. A semicolon will 
separate the upset from the conditioning event. P{U ; ei) 
denotes the probability of the upset U given that the edge 
ei is open, while P{U ; el) denotes the probability of the 
upset U given that the edge ei is closed. 

Using this notation and continuing the discussion pre- 
ceding the notational description, Pp[ABC;ei] ~ 0. 

Let each bond in a lattice be open with probability p, 
< p < 1. Let ei be an edge in E{T), and therefore in 
E{G), and let 62, 63, . . . , e„ be the other bonds in G. 
Let {Xj}^^-^ be independent uniform random variables 
on the interval (0,1). For each realization x of the 
{XjYj^^^, if Xj < p, call Cj open. Otherwise, ej is closed. 
Call this scenerio Model 1. Note that P^[ABC] denotes 
the probability that A, B, and C are all connected in G 
in this case. 

Consider the same lattice with the only difference 
being that ei is open with probability p + e. Call this 
scenerio Model 2, and note that Pp+g,p[ABC] denotes 
the probability that A, B, and C are all connected in 
this case. The other partition probability functions are 
defined in this case using the same subscript. 



Call Model 3 the same as Model 1 with p replaced by 
p + e in both the description and the notation. We use 
the same realization 'x in all three cases. 

Notice that if an edge is open in Model 1, it is also 
open in Model 2. So, if A, B, and C are connected 
through open bonds in Model 1, they are necessarily 
connected through open bonds in Model 2. Also, if 
an edge is open in Model 2, it is open in Model 3 and 
similar conclusions can be made. 

In what follows, P'^[ABC] is the probability that 
A, B, and C are connected in G only through T, 
that is, they arc not connected if any edge of T is 
closed. P'^IABG] is the probability that A, B, and 
C are connected in G — e, the graph G with edge e 
deleted, for some edge e e E{T). We then have that 
P^[ABC\ = P'^[ABC\ + P'^[ABC]. We obtain a lower 
bound on the derivative by considering the difference 
quotients: 

lim P?a^'^C]-P^[ABC] 
> lim ^°+e..[^ggl-^°[^ggl 

- £^0 

^ lim Pv+..MBC]+P^^^JABC]-(P^[ABC]+P^[ABC]) 
£^0 « 

Substituting 

Pj,+e,p[ABC] = {p+e)Pp[ABC;ei]+{l-p-e)Pj,[ABC;el] 



and Pp[ABC] = pPp[ABG;ei\ + {I - p)Pp[ABG;-ei\, 
we see that the last expression above equals: 
P^[ABC; e^]-P^[ABC;ei\+P^[ABC; er]-P^[ABC;el] 

Since P^[ABC; d] = p™"!, P^[ABC;eT] ^ 0, and the 
event that A, B, and C are connected in N with ei closed 
is contained in the event that the three are connected, 

P9,, [ABC] - [ABC] 
lim ^ > > 

c^o e ~ 

Similar reasoning is valid for the other partition 
probability functions, one of which will be given in 

detail here. Let P[AB] = P[ABC] + P[AB\Cl the 
probability that A and B arc connected in the lattice. 
Also, let ei be an edge in T on the unique path from A 
to B and let a be the number of edges in the path. Then, 

lim P°+e[^gl-P°[^gl 

£^0 « 

> lim ^p°+e,.[^gl-P°[^gl 

- £^0 « 

£^0 " 

^ lim <Pl[AB-M+P?[AB-M+Pj[AB-;e^]+P^[AB-;eI]) 

= 7f[AB; ei] + (P/ [AB; ei] + P/ [AB; el]) 
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> P' 



,a-l 



> 0, 



In this way, it follows that all upset probability func- 
tions for a generator of a lattice in the class have a pos- 
itive derivative in the interval (0,1). The exact same 
argument holds for the duals of these lattices, replacing 
p by q. Thus, the coefficient of e used in this method to 
determine equality of the bond percolation critical expo- 
nents is well-defined and positive. As a consequence, any 
lattice and dual lattice in the class described have equal 
values for the two critical exponents (3 and 7. 



V. EXTENSIONS AND GENERALIZATIONS 
A. Other Critical Exponents 

Additional percolation functions and critical exponents 
are considered in the literature. Define the two-point 
connectivity function by: 

T{p,x,y) = Pp[x y], where "a; ^ y" denotes the event 
that the sites x and y are connected by an open path. 
Define the correlation length by 



a?) 



oo\ 



where \y\ = max{|y(z)| ■.l<i<d},y= (y(l), . . .,y{d)). 

The critical exponents v, S, 77, and A^, for fc > 2, are 
given by the power laws: 



gp[|Cr;|C|<oo] 
Bp[|C|'=-i;|C|<oo] 



^{p) ~ \P^ Pc\^'' for some > 0, 

\p-Pc\^'^'', fc > 2, for some Aj. > 0, 

Ppa 1^ 1^*1 < 00] w n^/^ , n — > 00, for some 6 > 0, 

and, for a d-dimensional periodic graph, 

t{pc,v,x) « \x\'^^'^^^, \x\ 00, for some 77 > 0. 

In two dimensions, the values for the critical expo- 
nents defined in this paper are believed to be /3 = 5/36, 
7 = 43/18, S = 91/5, and ^ 4/3 (See Grimmett [13 
and Hughes Kesten [2,] proved that for a class of two- 
dimensional periodic lattices (assuming the limits defin- 
ing the exponents d and v exist), that 



S+1 

and that for all fc > 2 



2iy ^6 



1 



<5-f-l 



A. 



2v- 



1 



The class of lattices considered by Kesten includes the 
lattices in the class identified by Wierman and Ziff. 



We have already shown that f3{L) = f3{D) and j{L) = 
7(13) for any lattice L in the given class. Furthermore, 
since P^^ and are equal, T{qc,x,y) — T{pc,x,y) 
for all boundary vertices x and y. Thus, T{qc,x,y) 
and t{pc, x, y) decrease at the same exponential rate, so 
ri{K) — r]{M). Using Kesten's formulas, equality of the 
other critical exponents is evident. That is, 5{L) = S{D), 
v{L) = v{D), and A(L) = /\{D). 

Thus, the set of critical exponents are equal for a bond 
model and its dual in the given class. For example, the 
results apply to the lattice pairs mentioned in Ziff and 
ScuUard 7:|. Since each generator can appear in an infi- 
nite collection of self-dual hypergraphs, we have infinitely 
many lattices in which the set of critical exponents are 
equal. 

Our results do not establish any numerical values for 
the critical exponents. However, remarkable progress in 
this direction was made by Smirnov and Werner [5j| , who 
combined Kesten's scaling relations, knowledge of crit- 
ical exponents associated with the stochastic Loewner 
evolution process, and Smirnov's proof of conformal in- 
variance to determine the existence and values of critical 
exponents for the site percolation model on the triangular 
lattice. 



B. The Bowtie Lattice and Its Dual 

By splitting each vertical bond in the bowtie lattice 
into two bonds, each having probability 1 — ^1 — p of 
being open, Wierman [3] used the substitution method 
to determine the critical probabilities of the bowtie lat- 
tice and its dual. The method described in this paper 
was then applied to this dual pair, showing the equality 
of their bond percolation critical exponents. Using this 
idea, the method described in this paper is applicable to 
many other lattices. In fact, as long as a lattice is self- 
dual under the triangle-triangle transformation (see Ziff 
and ScuUard Q, and Wierman and Ziff [9^]), functions 
of p can be assigned as edge probabilities of a given lat- 
tice. That is, we can assign to each edge of the generator 
the probability fe(p) of being open, where fe(p) is an 
increasing, right-continuous function of p. Under these 
conditions, repeating the main argument of the paper 
gives that, for any edge e in the minimal tree connecting 
A, B, and C, 



lim 



P^_,,[ABC]-P.^[ABC] 



Pi^[ASC|e]-Pi^[Ai3C|e]) 



and it is clear that the coefficient of e in the expression 
for (5(e) is well-defined and positive. Equality of bond 
percolation critical exponents for such lattices and their 
duals follows from this fact. 
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C. Site Percolation 

Since the bond problem of a lattice is e quiv alent to the 
site problem on its line lattice [1, the result 

described in this paper applies to the collection of site 
problems obtained by reformulating the bond problems 
as site problems. As a result, we have a collection of site 
problems in which any lattice belonging to the collection 
shares the same value of its critical exponents with its 
matching lattice. So, for the bond model on any given 
lattice in the class identified by Wierman and Ziff, we can 
immediately identify three other lattices with the same 
values for the set of critical exponents: the bond model 
on the lattice's dual, the corresponding site problem on 
the line lattice, and the site problem on the line lattice 
of the dual. Since each generator can appear in an infi- 
nite collection of self-dual hypergraphs, we have infinitely 
many site problems in which the set of critical exponents 
are equal. _ 

The martini. A, and B lattices discussed in ScuUard 
are not line graphs of underlying lattices, so the results 
in the previous paragraph do not apply to the site per- 
colation models on these lattices. It is plausible that the 
approach of this article may be applicable to such site 
percolation models, but that has not yet been shown to 
be valid, and is a subject of further study. 

VI. CONCLUDING REMARKS 

The method used by Wierman to show equality of the 
bond percolation critical exponents for the triangle and 



hexagonal lattices was used to show the equality of the 
bond percolation critical exponents for the martini lat- 
tice and its dual. This computational proof was then 
extended to show that a lattice constructed from an in- 
finite connected planar periodic 3-uniform hypergraph 
with one axis of symmetry, using a generator which is 
a finite connected planar graph with three boundary ver- 
tices, has the same values for its bond critical exponents 
as its dual, thereby generalizing the results of Wierman 
[Sj. Since the mentioned class of lattices is infinite, there 
are infinitely many lattices that have the same bond per- 
colation critical exponents as their duals. Moreover, since 
using the same generator on different self-dual hyper- 
graphs does not affect the computations, there are an 
infinite number of bond and site models for which the 
critical exponents are the same. This result gives mathe- 
matical evidence that the values of the critical exponents 
may only depend on the dimension of the lattice, sup- 
porting the universality hypothesis. Note that the re- 
sult does not say that the bond percolation critical expo- 
nents have the same value for all lattices in the mentioned 
class. Using different generators produces different upset 
probability functions, so no relations for the bond per- 
colation critical exponents between generators has been 
determined. If two non-isomorphic, non-dual generators 
were discovered that had equal upset probability func- 
tions at criticality, the result could likely be extended to 
give equality of bond percolation critical exponents for 
these lattices. The authors of this paper have thus far 
been unable to identify two such generators. 
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